Polaron contributions to the biexciton binding energies in self-assembled quantum 

dots 



(N 

o 

(N 

a 

ON 



c 

o 

(N 
> 
O 
VO 
oo 
cn 

On 
O 
O 



X 



Pawel MachnikowskQ 

Institute of Physics, Wroclaw University of Technology, 50-370 Wroclaw, Poland 

The contribution to the biexciton binding energy in quantum dots resulting from the interaction 
with longitudinal optical phonons is estimated by performing the configuration-interaction calcula- 
tion of the few-particle states in a simple model of the confining potential and including the phonon 
corrections by means of a perturbation theory. It is found that the polaron contribution tends to 
compensate the Coulomb-related biexciton shift (binding energy) and reduces its value by several 
to even 30%, depending on the material parameters of the system. 
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The biexciton binding energy (also referred to as the 
biexciton shift) is one of the most important parame- 
ters determining the optical properties of quantum dots 
(QDs). For instance, the spectral shift between the ex- 
citon and biexciton transition makes it possible to inde- 
pendently address various degrees of freedom of the biex- 
citon system which opens the way to a range of quantum 
optical control schemes^. The magnitude of this shift 
is usually of the order of a few meV while its sign may 
be either positive or negative 3,4 . Theoretically, the ef- 
fects of Coulomb correlations in few-electron QDs, un- 
derlying the biexciton shift, can be calculated using a 
configuration-interaction scheme^ either using a simple 
model of confinement^— or built upon more realistic wave 
functions obtained within the multiple-band kp or pseu- 
dopotential theory 9 . 

One factor that is missing from all the existing cal- 
culations is the effect of the interaction between car- 
riers and longitudinal optical (LO) phonon which is 
of particular importance in confined system o 10 ' 11 and 
leads, among other effects, to polaron shifts and res- 
onant anticrossings 1 ^— . According to theoretical cal- 
culations, the magnitude of these features may be as 
large as a few meV for a single confined exciton in a 
self- assembled Q D 15 i 16 . Since this is comparable to the 
biexciton binding energy it seems interesting to study the 
phonon-related shift for a confined biexciton. 

The purpose of this paper is to estimate the polaron 
contribution to the biexciton binding energy in a QD. 
This will be done by calculating the polaron correction 
to the exciton and biexciton states using the perturbation 
theory. The few-particle states themselves will be calcu- 
lated within the usual configuration-interaction scheme 
with a simple parabolic model of confinement. 

The quantum dot will be modeled by axially symmet- 
ric harmonic oscillator potentials for electrons and holes 
with a strong confinement in the growth direction. In 
order to mimic (in a simple way) the spatial charge sep- 
aration observed in many QD structures, the minima of 
the potentials for electrons and holes are allowed to be 
displaced along the growth direction and located at dif- 
ferent points z e and Zh, respectively. The single particle 



Hamiltonian is then 
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where "e/h" corresponds to the electron and hole wave 
functions, m / h are the effective masses, and w e /h) ^o/h 
are the harmonic oscillator frequencies for the confine- 
ment potential in the xy plane and along the z axis, re- 
spectively, with £> e /h 3> w e /h- The motion along the z 
axis (growth direction) will be restricted to the ground 
state. It will be assumed that the confinement width 
of a single (non-interacting) particle is the same for the 
electrons and holes, 
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for the in-plane and axial confinement, respectively. 

The calculations will be performed in the basis of 
single-particle wave functions 
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is the wave function in the z direction and 
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is the Fock-Darwin wave function of the 2-dimensional 
harmonic oscillator—, where C'™ '(x) is the Laguerre 
polynomial. Here n = 0,1,... is the radial quantum 
number and m = 0,±1,... is the angular momentum 
quantum number. The length parameters Z c and lh are 
determined by minimizing the variational exciton ground 
state 
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where 



exciton energies and the corresponding eigenstates 
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eo is the vacuum permittivity, and s s is the static relative 
dielectric constant. In the following, we will use a single 
Greek index a, /3, . . . to represent the pair of quantum 
numbers (n, to). 

The confined few-particle electron-hole subsystem is 
then described by the Hamiltonian 

after after 

+ Y Y Vcc ( a 'P ;a '>P^ a ''°'<r a l°' a f :i '< 7 ' aa ' a 

CK0CK 1 'ft' <7<7' 

afta' 'ft' crcr' 

- E E^( a ># a ^>W>^' a «'<#) 

where a aCT , a* acr and ft, acr , /i^ CT are the annihilation and cre- 
ation operators for electrons and holes, respectively, cor- 
responding to the states described by the wave functions 
given by Eq. (jlj and spin orientation a, 

and 
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with i, j = e, h. 

The interaction between the carriers and LO phonons 
is described by the Frohlich Hamiltonian 16 
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where e = (l/e„— l/£oo) _1 is the effective dielectric con- 
stant (eoc is the high-frequency relative dielectric con- 
stant), is the frequency of the LO phonons (assumed 
constant in the narrow range of the relevant wave vec- 
tors), bk and are the annihilation and creation opera- 
tors for the LO phonon with the wave vector k, V is the 
normalization volume for phonons, and the form factors 
are given by 

J^\k) = J d 3 r^ h >(ry k - r * { ; /h) (r). 

First, the Coulomb part of the biexciton shift AEq is 
determined. To this end, the Hamiltonian ([2} is diagonal- 
ized in the truncated basis of one-pair states yielding the 
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where |0) is the empty dot state. The subtle fine struc- 
ture effects are not included in the model, so the states 
are degenerate with respect to spin configurations. As 
both the Coulomb interaction and the carrier-phonon 
couping conserves separately the electron and hole spins 
and we are interested in the ground state only it is suf- 
ficient to take into account the two-pair states with sin- 
glet spin configurations. The energies of the biexciton 
states are therefore found by diagonalizing the Hamil- 
tonian -ffe-h in the basis of two-pair spin-singlet states, 

(2) 

from which one finds the energies E„ and the corre- 
sponding eigenstates 
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where j" and J, denote spin orientations, r) a p = (1 — 
\/2)<5 Q/ 3 + y/2, and comparing the labels a, j3 is performed 
in the sense of an arbitrary but fixed ordering of the 
single-particle states. The Coulomb part of the ground 
state biexciton shift is then simply 
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Next, one has to compute the phonon- induced correc- 
tions to the exciton and biexciton states. Even though 
the polaron shifts for single carriers can be quite large, 
those appearing for globally neutral exciton states are 
smaller due to partial cancellation of the electron-phonon 
and hole-phonon interactions. Therefore, a reliable esti- 
mate of the polaron shift in the cases considered here can 
be obtained by means of the second order perturbation 
theory. The calculation is done again separately for the 
exciton and biexciton cases. 

In the case of an exciton, one rewrites the carrier- 
phonon interaction in terms of the one-pair eigenstates 
[Eq. (for an arbitrary, fixed spin configuration), 
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In an analogous manner, for the biexciton, one writes 
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TABLE I. Material parameters used in the calculations, mo 
is the free electron mass. 
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The phonon-induced correction (polaron shift) to the 
ground state energy is then given by 
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where z = 1,2 refers to the exciton and biexciton cases, 
respectively. The phonon correction to the biexciton shift 
is 



AE ph = 6E<® - 26E ( - 1 \ 
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The total biexciton shift is AE = AE C + AE ph . 

In the numerical calculations, two sets of parameters 
will be used, representing two groups of material sys- 
tems: GaAs (a moderately polar III-V compound) and 
CdTe (a II- VI compound with a more polar character) 
(see Tab.[T]). The numerical computations are performed 
using the basis of 6 electron and hole shells (2|m|+n < 5), 
that is, 21 electron and hole wave functions. This is 
sufficient to assure the convergence of numerical results 
within ~ 0.1 meV of possible error, which is sufficient in 
view of the much larger magnitude of the discussed effect 
(a few meV). 

The results for the first set of parameters are shown 
in Fig. [I] where the dependence of the biexciton shift on 
the e-h displacement D = z c ~ z n is shown for two lateral 
QD sizes. Although the polaron contribution does not es- 
sentially change the qualitative behavior of the biexciton 
shift it constitutes an important correction to the over- 
all value. For this material system, the phonon-induced 
contribution amounts to roughly 15-17% of the Coulomb 



FIG. 1. The biexciton shift (black solid lines) for a GaAs- 
like QD system, as well as its Coulomb (red dashed lines) and 
phonon (blue dotted lines) components for the two lateral QD 
sizes characterized by hw B — 40 meV (a) and 70 meV (b). In 
both cases l z = 0.2Zn- 
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FIG. 2. As in Fig. [T]and for the same two values of hw e but 
for a CdTe-like system. The gray dash-dotted line shows the 
biexciton shift with the correction due to the axial shift. 



part over the whole range of e-h displacements D stud- 
ied here, except for the region where the binding energy 
changes sign and this ratio is poorly defined. The rela- 
tive value of the polaron part is very similar for both QD 
sizes even though the absolute magnitude of the total 
binding energy is larger for the smaller dot (larger hio e ). 
It turns out that the phonon correction has the opposite 
sign to the Coulomb contribution and, therefore, reduces 
the absolute value of the binding energy. Interestingly, 
both these contributions change sign at almost the same 
value of D. 

In Fig. [51 analogous results for a more strongly polar 
system are shown. The magnitude of both the contri- 
butions to the binding energy is now much larger. Also 
the role of the phonon correction becomes more impor- 
tant and its value reaches almost 30% for both dot sizes. 
Again, the signs of the two contributions are opposite 
over most of the parameter range studied. In this case, 
the sign change still appears in the same region for both 
contributions but not as closely as in the previous case. 

In the presented discussion, the wave functions along 
the growth direction were restricted to the ground state, 
that is, the charge distribution was supposed to be rigid 
with respect to axial shifts. However, for D > 0, the 
Coulomb attraction will lead to an axial shift of the elec- 
tron and hole charge distributions, which will contribute 
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to the biexciton shift. Accounting for this effect requires 
extending the computational basis by including the first 
excited state wave functions along z. By virtue of the 
Ritz theorem, the same goal (to the leading order) can 
be achieved by a variational minimization of the system 
energy with respect to the positions of the centers of the 
wave functions along z for a given distance D between 
the minima of the electron and hole confinement poten- 
tials. The gray dash-dotted line in Fig. [2] shows the biex- 
citon shift for the CdTe-like system including the addi- 
tional contribution from the axial shift obtained by such 
a procedure (with fixed l e i^ found earlier). The correc- 
tion reaches almost 1 meV at D = 2 nm, which is no- 
ticeable, although not particularly large compared to the 
values without the correction. For a GaAs-like system, 
the correction does not exceed 0.2 meV and is therefore 
relatively much smaller. The shift of the charge distribu- 
tions does not exceed 0.04 nm and 0.1 nm for the GaAs 
and CdTe system, respectively. Therefore, the correction 
to the phonon part, which is only sensitive to the ge- 
ometry of the wave functions, will reach 0.03 meV and 
0.25 meV, respectively, as can be deduced from the slopes 
of the phonon curves in Figs. [T] and [2] In both cases, this 
corrections is rather small compared to the overall values 
of the biexciton shift. 

Finally, it can be argued that the results presented 
above are in a way nontrivial, at least in the sense that 
they cannot be reproduced using a simple model. In- 
deed, if one considers only exciton and biexciton wave 
functions in a product (Hartree) form and takes into ac- 
count that the wave functions for the two electrons (or 
for the two holes) that differ only by their spin orienta- 
tion must be the same then the resulting biexciton shift 
must be positive even if the single-particle wave func- 
tions in the biexciton state are allowed to be different 
from those in the exciton state. Therefore, the negative 



value of the shift (positive binding energy) is a purely 
correlation-induced effect. 

For the phonon effect, it seems impossible to draw such 
a general conclusion since the polaron correction depends 
strongly on the charge cancellation and hence on the ex- 
act shape of the wave functions which can differ between 
the exciton and biexciton cases. However, if one restricts 
the discussion to a single-orbital model in which only one 
wave function is available for electrons and one for holes 
(the same in the exciton and biexciton configurations) 
then it can be seen from Eqs. © and (|7|) that 

g (2) (k) = 2 \T {h \k) - T ic) (k)\ = 2g (1 \k), 

where the indices of the form factors have been omitted. 
As there is only one state, only one term (n = 0) remains 
in Eq. © and, clearly, W? (2) = 4<5i? (1) . Therefore, ac- 
cording to Eq. ©, AE'ph = 26E& < 0, since the ground 
state polaron shift is negative. Clearly, this result is op- 
posite to that obtained from the full model, where the 
phonon correction to the biexciton shift is positive for 
small D. Therefore, exact shape of the wave functions, 
including admixtures of higher single particle states is 
essential for correct modeling. 

In summary, the LO phonon-related (polaronic) contri- 
bution to the biexciton shift (biexciton binding energy) 
in quantum dots has been calculated. It turns out that 
this correction is about 15% of the dominating Coulomb 
contribution in a typical III-V material and reaches al- 
most 30% of the Coulomb term for a more strongly polar 
II- VI system. The phonon correction usually has an op- 
posite sign to the Coulomb part and, hence, reduces the 
binding energy of a confined biexciton. 
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